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Abstract 



We study a system of nonlinear partial differential equations resulting 
from the traditional modelling of oil engineering within the framework of 
the mechanics of a continuous medium. Existence and regularity of the 
optimal solutions for this system is established. 
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1 Introduction 

We are interested to the existence and regularity of optimal solution for the 
following "dead oil isotherm" problem: 



where Q is an open bounded domain in R 2 with a sufficiently smooth boundary. 

Equations l|l.lf> serve as a model of an incompressible biphasic flow in a 
porous medium, with applications in the industry of exploitation of hydrocar- 
bons. To understand the optimal control problem we consider here, some words 
about the recovery of hydrocarbons are in order. At the time of the first run of 
a layer, the flow of the crude oil towards the surface is due to the energy stored 
in the gases under pressure or in the natural hydraulic system. To mitigate 
the consecutive decline of production and the decomposition of the site, water 
injections are carried out, before the normal exhaustion of the layer. The water 
is injected through wells with high pressure, by pumps specially drilled with 
this end. The pumps allow the displacement of the crude oil towards the wells 
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d t u — Aip(u) = div (<?(u)Vp) 



inQ T = nx (0,T), 
in Q T = nx (0,T), 



d t p - div (d(u)Vp) = / 
P\dn = °> Pl i= o=Po, 



(1.1) 



of production. The wells must be judiciously distributed, which gives rise to a 
difficult problem of optimal control: how to choose the best installation sites 
of the production wells? The cost functional to be minimized comprises all the 
important parameters that intervene in the processes. 

Existence and uniqueness to the system , for the case when the term dtp 
is missing but for more general boundary conditions, is established in Opti- 
mal control of systems governed by partial differential equations is investigated 
in literature by many authors, we can refer to 0E1II5]. To study existence and 
regularity of solutions which provide Gateaux differentiability of the nonlinear 
operator corresponding to we are forced to assume more regularity on 

the control / as well as to impose compatibility conditions between initial and 
boundary conditions. The considered cost functional comprises four terms and 
has the form 

J(u,p, f) = \ \\u - u\\1 Qt + ~ |b - p\\1, Qt + § \\f\\ll. QT + y WWWIqt 

(1.2) 

where 1 < qo < 2, j3\ > and 02 > are two coefficients of penalization; 
U and P are given data. Here u is the reduced saturation of the phase oil 
at the moment t. The initial saturation is known and p is the total pressure. 
The first two terms in 1|1.2|) make possible to minimize the difference between 
the reduced saturation of oil and a given U , respectively the global pressure 
and a known initial pressure P. We remark that the choice of the objective 
functional is not unique. We can always add further terms of penalization to 
take into account other properties which one may want to control. The paper is 
organized as follows. In Section El we set up the notation, the functional spaces 
and some important lemmas used throughout the work. Section [3] is devoted 
to the existence of optimal solutions. We obtain necessary estimates on the 
sequence minimizing the cost functional which allows us to pass to the limit. 
Finally, in Section 01 we establish a regularity theorem. 

2 Notation and Functional Spaces 

In the sequel we suppose that tp, g and d are real valued C^-functions satisfying: 

(HI) < ci < d(r), tp{r) < c 2 ; \d'(r)\, \if'{r)\, \<p"(r)\ < c 3 Vr e R. 

(H2) u , po £ C 2 (Q), U, P £ L 2 {Q T ), where u , p : Q -> K, U, P : Q T -> R, 
and u \ dn = p \ aQ = 0. 

We consider the following spaces: 

W^°(Q T ) := LP (0,T,Wfa)) ={ue L P (Q T ), Vu e 27(Q T )} , 
endowed with the norm ||u|| w i.o ( q t) = |M| P:Qr + ||Vw|| PiQr ; 

W^iQr) := {u e W^°(Q T ), V 2 u,8 t u e L^{Q T )} , 



2 



with the norm IMI^.i^ = IMI^.o^ + ||V 2 u|| PjQt + \\d t u\\ p ^ T ; 



V := 



{u £ W^°(Q T ), d t u £ L 2 (0,T,W 2 -\n))} . 



We now state some important lemmas that are used later. Lemma 12.11 is 
needed in the proof of our existence result. 

Lemma 2.1 ([T]). A ssume n C R ra is a bounded domain with a C 1 -boundary, 
and a matrix A(x,t) = (Aij(x,t)) satisfying the conditions 

3 7o > such that Aij (x, > 7o |£| 2 V£ £ M™ , 

Aij £ L°°(Q T ) , ,4^=%. ( ' J 

Assume also f £ L 2<Io (Qt), uq £ WagoO^O / or some <Zo > 1 arl ^ u £ 
C ([0, T]; L 2 (f2)) fl Wq' (Qt) ^ e a weaA; solution to the equation 

dtu — div (A(x,t)\Iu) = f in Qt , 

"1(50 = °> M U=0 = "0 ■ 

Then, there exists a constant q > 1, depending on n, qo, jq, Qt, and H-AH^ q t , 
such that u £ W^ q (Qt) and the estimate 

l|V«|| a?i Q T <c([[/|| 25igT + ||« [| w -i <( n)) 

holds. 

We use the following two lemmas to get some regularity of weak solutions. 

Lemma 2.2 (De Giorgi-Nash-Ladyzhenskaya-Uraltseva theorem [6 ). 

Assume Qt = f2x (0, T), O C R" a C 1 -bounded domain, and let f £ L s ' t {Qt) — 
L s (0,T,L r {Q)), u £ C a (Cl) for some a Q > 0, u \ gn = cmd 
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r + 2^ <1 ' 

Assume (|2.1|l Zio/ds and Zet u £ W^' (Qr) &e a weafc solution of l|2.2|l . Then, 
there exists a > swc/i /j/iai u £ C a ^ {Qt) and 

IHc"-* (Q r ) ^ C (II/Hl-XQt) + H U 0|lc°(O)) • 

Lemma 2.3 (0). For any function u £ C a -t (Q T )ni 2 fo, T; M^ 1 (O) n Wf(fi) 

there exist numbers No, go such that for any g < go there is a finite covering of 
by sets of the type fl g (xi), Xi £ Q, such that the total number of intersections 
of different &,2g(xi) = fin B2 g (xi) does not increase Nq. Hence, we have the 
estimate 

\\Vu\\i QT < c \\uf ca , nQT) e *° (||v 2 u||^ r + 1 HVulli^) . 
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3 Existence of Optimal Solution 

We denote by (P) the problem of minimizing l|1.2f> subject to in the class 
(«,P,/)e^ 1 (Qr)xyxL 2 (Q T ). 

Theorem 3.1. Under hypotheses (H1)-(H2) there is a q > 1, depending on 
the data of the problem, such that there exists an optimal solution (u,p, fj of 
problem (P) verifying: 

ueW 2 > x {Q T ), 

peC([0,T];L 2 (Q))nW^ (Q T ), d t p £ L 2 (0, T, W 2 _1 (fi)) , 
/6£ 2 «(Q T ), d t feL 2 (Q T ). 

Proof. Let (u m ,p m , f rn ) € W 1 ' (Qt) x V x L 2q °{Q T ) be a sequence minimizing 
J(u,p, /). Then we have 

(f m ) is bounded in L 2qn (Q T ), 

(d t f m ) is bounded in L 2 (Q T ). 

Using the parabolic equation governed by the global pressure p and Lemma I/TT1 
we know that there exists a number q > 1 such that 

iiv P ™ii 2g , Qr < (um 2g,Qr + ll u || W\ q (fi) ) ■ 

Multiplying the second equation of by p, using the hypotheses and Young's 
inequality, we get 

sup|b m ||L, + ||Vp m ||^ QT < C ||/ m ||^ QT . 

Furthermore, we have that d t p m is bounded in L 2 (0,T;W 2 1 {n)). By Aubin's 
Lemma 8 , (p m ) is compact in L 2 (Qt)- Using now the first equation of (|l.lfl 
we have 

d t u m - (p'(u m )Au m - ip"{u m )\Vu m \ 2 = div{g{u m )Vp m ). 

Hence 

\\u m \\ w 2, 1{QT) <c, 

where all the constants c are independent of m. Using the Lebesgue theorem 
and the compacity arguments of J. L. Lions 8 we can extract subsequences, 
still denoted by (p m ), (u m ) and (f m ), such that 

p m ^p weakly in W£ (Q T ), 

d t p m -> <9 t p weakly in i 2 (0,T; Wf^Q), 
p m — > p strongly in L 2 {Qt), 
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p m ^pa.e. in L 2 (Q T ), 

u m a.e. in L 2 (Q T ), 

f m ^ / weakly in L 2 «°(Q T ), 

dtf m -> 3*7 weakly in L 2 (Qt). 

The existence of an optimal solution (u,p, f) follows, in a standard way, by 
passing to the limit in problem and by using the fact that J is lower 

semicontinuous with respect to the weak convergence. □ 

4 Regularity of Solutions 

We now prove some regularity to the solutions predicted by Theorem l3.ll 

Theorem 4.1. Suppose that (HI) and (H2) are satisfied and let (u,p, f) be an 
optimal solution of our problem (P). Then, there exist a a > such that the 



following regularity conditions are verified: 

p e C a <% (Q T ) , (4.1) 

u,pe Wl fi {Q T ) , (4.2) 

u,peW^(Q T ), (4.3) 

d t u, d t peL™ (0,T;L 2 (n))nW}'°(Q T ), (4.4) 

ueC*(Q T ), (4.5) 

ueW 2 +(Q T ), peW 2 q l(Q T ). (4.6) 



Proof. First, we remark that 1|4.1[) is an immediate consequence of Lemma \'2. 21 
To show the other results, we begin by proving the following lemma. 

Lemma 4.2. Consider (u,p,f) solution of (11.1(1 . Assume that hypotheses (HI) 
and (H2) hold. Then, 

sup ||Vp|£ n + || V 2 p||' < c (\\V P \\% Qt + || Vu\\1 Qt ) + c 
te(o,T) x 

where c depend on uq and f . 

Proof. From the second equation of we have 

d t p - d(u)Ap = d'{u)S7u.S7p + f . 
Multiplying this equation by dtp and integrating over Q, we obtain 

\\9tp\\l + ~ || Vp|| 2 < c / \WpVud tP \ dx+ f \fd tP \ dx . 

Using Young's inequality and integrating in time, we get the desired estimate. 

□ 
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To continue the proof of Theorem 14.11 we need to estimate ||Vw||4,q t in 

function of ||Vp||4,g r . Then, taking into account the first equation of it 

i i 

is well known that u € W A ' 2 (Qt) and 

||V«|| 4 ,q t < c\\\Jp\U, QT (4.7) 
(see 0]). Using Lemma \'2. 'M we have that for any g < go 

\\V P \\i QT < c bl|*„. f {Qt) g 2a {||V< Qt + 1 ||Vp||^ T } + C„ 0iPO)/o . 

Calling Lcmma l2~2l we then get (14.21) for an eligible choice of g. After using (|4.7|l 
we obtain that u G W4' (Qt)- On the other hand, we have by the first equation 
of and l|4.2jl that u E W^iQr)- Moreover, it follows by Lemma FQ1 and 
the fact that u e W^ l {Q T ) that p 6 W^'^Qt). 

Now, in order to prove (|4.4|) . we differentiate both equations of with 
respect to time: 

dttu — div ((p' '(w)V dtu) — div (tp" (u)V dtuV u) 

= div (g'(u)d t uVp) + div (g(u)Vd t p) , (4.8) 

dup - div {d{u)Vd t p) - div {d! (u)Vd t uVp) = d t f. (4.9) 
Multiplying (|4.9(l by d t p and integrating over Q we get 

^liatP||a,n + cliatVpH^n < cj + c||%>||i,n + c ^ |d tU V P V9 tP | dx. 
By Young's inequality we have 

/ |%tVpV%>| da; < ||ft«Vj»|| 2 ,n||ftVp||a,Q 

<c||9 t uVp||| n + |]|9 t Vp||| in . 
On the other hand, by Holder's inequality we obtain 

\\d t uVp\\l n = [ |cH 2 |V?f 




= \\dtu\\l n \\Vp\\l 



Using the following multiplicative inequality jS] 

\\dtu\\ln < c|[ft«||2,n||ftVu|| 2 ,n Vu e W^fi), 
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we obtain 



J \d t uV P Vd tP \ dx < c||^tt||2 )n ||Vp||t,n + |l|9tVp||| Q 



< c||au|| 2|n ||ftVu|M|Vp||i, n + l\\d t Vv\\l,n 

< cllfttill^nllVpHtn + o\\d t Vu\\l n + |||^Vp||l^. 

Then, 

§- t \\d t p\\in + 4dtVp\\ln 

<c f + c\\d tP \\l n + c\\d t Vu\\l n + c\\d t u\\l n \\\7p\\l n . (4.10) 
Multiplying (|4.8(l by d t u and integrating over ft, we get 

§- t \\dtug >n + c\\d t Vu\\l n 

< - ip"(u)d t uVud t Vu - / g'(u)d t uVp\/u - / g(u)d t Vp\/u. 
Jn Jn Jn 

Similar as before, we have: 

1 <c\\d t u\\ij\Vu\\jn + c\\d t Vu\\l n , 



ip"(u)d t uWud t Vu 
g'(u)d t u\7p\7u 



<c\\d t u\\ln\\Vp\\jn + c\\Vu" 2 



2,f2i 



g{u)d t V P Vu < c\\d t Vp\\l n + c\\Vu\\l 

Jn 

It follows by using l|4.7[l that 

§- t \\dtu\\l n + c\\dtVu\\l n < c||ft«||i in ||Vp||J )n +c||ftVi>||i in +c||V«||i in . (4.11) 
Calling l|4.10jl and l|4.11|l together, it yields 



j t {\\dtu\\ln + \\d tP \\l n } + \\dtVu\\l n + \\dtVp\\l n 

We thus obtain (|4.4I) by applying Gronwall lemma. On the other hand, we have 
||<9*M||4,n < c\\dtu\\2,n\\dtVu\\ 2: n, 

and from l|4.4|l we obtain <9 t u 6 L a (Qt)- Using then 14.2fl and the fact that 
WI(Qt) c — > C*(Qt), the regularity estimate l|4.5|l follows. Finally, the right 
hand side of the first equation of 11.1(1 belongs to L a (Qt) L 2q °(QT) as 
2q < 4. Using thus |@3J| we get u € W$£(Qt). Since / e L 2 9°(Q T ), the 
same estimate follows for p from the second equation of the system (|1.1|) and 
we conclude with l|4.6|> . □ 
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